A phantom force induced by the tunnehng current, characterized on Si(lll) 
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Simultaneous measurements of tunneling currents and atomic forces on surfaces and adsorbates 
provide new insights into the electronic and structural properties of matter on the atomic scale. We 
report on experimental observations and calculations of a strong impact the tunneling current can 
have on the measured force, which arises when the resistivity of the sample cannot be neglected. 
We present a study on Si(lll)-7x7 with various doping levels, but this effect is expected to occur 
on other low-conductance samples like adsorbed molecules, and is likely to strongly affect Kelvin 
probe measurements on the atomic scale. 



Scanning tunneling microscopy (STM) sparked enthu- 
siasm in scanning probe microscopy with images of the 
adatoms of Si(lll)-7x7 [1 . The atomic force micro- 
scope (AFM) removed the requirement for a conducting 
substrate [2] but brought more than just the possibility to 
measure on insulators [3]: Subatomic [4] and submolec- 
ular [5 imaging have also been demonstrated. These 
successes have brought strong interest in combined STM 
and AFM (e.g. [6j), however the independence of force 
and current measurements remains an open issue [THTT]. 

Frequency- modulation AFM (FM-AFM) is a technique 
in which the interaction between tip and sample is mea- 
sured by the frequency shift, A/, of an oscillating tip 
from its eigenfrequency, /o [12 . A/ can be formulated 
as a measure of the force gradient, kts = — where z 
is the distance from the surface. A/ is also a function of 
the spring constant of the oscillator, /c, the amplitude of 
oscillation, A, and and can be approximated at small 
amplitudes by A/ « {fo/2k)kts [13]. In short, a de- 
crease in A/ indicates that the force between the tip and 
sample is becoming more attractive. 

In contrast to the tunneling current, /, measured with 
STM, A/ is not monotonic as a function of z. The lo- 
cal tip-sample interaction is usually well-represented by 
a Morse potential from which the force and the force gra- 
dient can be derived, as shown in Fig. [TJa). In region I, 
kts decreases as z decreases. Long-range forces (e.g. van 
der Waals) cause attractive interaction between tip and 
sample. It is in this region STM is usually conducted 
on Si(lll)-7x7, at setpoints under 10 nA at IV, corre- 
sponding to tip-sample distances greater than 5 A [14 . 
In region II, kts increases as z decreases. The waveform 
overlap between tip and sample causes measurable en- 
ergy increase due to Pauli repulsion, which states elec- 
trons may not occupy the same quantum state [15]. 

In this Letter, we report upon the effect of bias voltage 
on FM-AFM of Si(lll)-7x7. At tip-sample distances 
corresponding to normal STM setpoints, one expects a 
decrease in frequency shift as the tip moves laterally with- 
out feedback over an adatom, due to the increase in at- 
tractive force m]. However, with the application of a 
moderate bias voltage (>1.0V), one is able to observe a 
frequency shift increase as the tip moves over an adatom. 



Moreover, FM-AFM images taken with this applied bias 
voltage can show atomic contrast at tip-sample distances 
300 pm further from the surface than is required to image 
with no applied bias. We propose a model incorporating 
sample resistance where the observed frequency shift is 
caused by a decrease in the electrostatic attraction be- 
tween tip and sample. 

Experiments were performed with a qPlus sensor with 
k = 1800 Nm-^ Data were collected in constant height 
mode with both a home-built microscope operating at 
room temperature in UHV and, where explicitly stated, 
at 4.2 K with an Omicron LT-SPM. Two types of Si(lll) 
samples were used: a high-doped sample corresponding 
to a resistivity p = 0.010 — 0.012 Q cm at 300 K and a low- 
doped sample with p = 6 — 9Q cm at 300 K. Si(lll)-7x7 
was prepared with several flash and anneal cycles. 

Figure [TJb) to (e) show simultaneously acquired / and 
A/ data of the low-doped Si(lll) sample. In (b) and (c), 
the tip bias was Vnp = —1.5 V. STM data show the clear 
structure of the 7x7 reconstruction, with all adatoms in 
the unit cell having approximately the same intensity. 
The A/ data show an increase in frequency shift above 
adatoms. In (d) and (e), the tip bias was Vtip = +1.5 V. 
The STM data show different features now in the outlined 
unit cell: over adatoms in the faulted half (the lower six 
adatoms), we record greater absolute current than those 
in the unfaulted half, as expected [17 . The A/ data, 
however, also have stronger contrast above the faulted 
half unit cell. In a simple picture of AFM in which the 
total electron density is measured [18 , one would not 
expect this A/ contrast to depend upon bias voltage. 
Also, while one might initially propose that this increase 
in A/ over adatoms indicates that we are in region II of 
the Morse potential, we show later (with respect to the 
data in Fig. [3| that this is not the case. 

The relation between A/ and / channels can be fur- 
ther characterized. For each pixel in image Fig. [ijb) to 
(e), force and current data have been acquired. Before 
a pixel- by-pixel comparision of / and A/ can be made, 
however, the low bandwidth of the PLL must be taken 
into account. At scan speeds even as low as 20nms~^, 
it can cause a noticable lateral displacement between A/ 
and / data. Consider a scan line 10 nm long with 256 
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FIG. 1. a) From a Morse potential, the force and force gra- 
dient can be calculated. At Vtip = —1.5 V, simultaneous b) 
current and c) force data are collected without / or A/ feed- 
back. Similar data are collected at Vtip = +1.5 V; the d) 
current data appears to be inverted because in opposite bias 
voltages, current flow is reversed, however e) A/ again in- 
creases above adatoms. Data were collected at A = 400 pm, 
/o = 25 908 Hz. In STM data, a unit ceh of 7x 7 is highlighted; 
images are lOnmxlOnm. 



pixels. The bandwidth of our PLL is 120 Hz. Assuming 
the / data to be instantaneous, the A/ data will be offset 
by (256 pixels/ 10 nm) • 20nms~^ / 120 Hz ^ 4 pixels. In- 
dependently, a cross correlation of Fig. [TJb) and (c) show 
a 4 pixel offset. 

Each measurement of / thus has an associated A/ mea- 
surement, and these are plotted (for the A/ and / data 
shown in Fig. [l] (b) and (c)) in Fig.j^Ja). The fact that 
there is some correspondance does not come as a surprise, 
as regularly force and current images of the same surface 
appear similar. What is surprising is that the data are 
quite linear and yield an increasing frequency shift over 
the entire current range. 

The slope of the linear fit to the data in Fig. |2|a) 
is a measure of the A/ response as a function of /. 
This analysis was repeated with images at various bi- 
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FIG. 2. a) From Fig[l];b) and (c), a comparison between / 
and A/ can be made. Their relationship is linear; the line 
shown has a slope of 4.36Hz/nA. b) Following similar anal- 
yses at other biases (Vtip=-1.5, -1.0, 1.0 and 1.5V) the slopes 
of the fits can be plotted as a function of the bias. / has the 
opposite sign as Vtip. A linear fit is shown as a guide to the 
eye. 



ases (Vtip = -1.5 V, -1.0 V, 1.0 V and 1.5 V). Fig. 3(b) 
shows the slopes of these fits with respect to applied bi- 
ases. The observation is that not only does this repulsive 
A/ signal scale linearly with current, but that this A/// 
slope itself scales with applied bias voltage. 

We then performed measurements with the high-doped 
sample. At room temperature, this A/ contrast was not 
observable in the current range of |/| < 5nA. At 4.2 K, 
however, it was, as shown in Fig. [sj^a): In cooling the 
system, the effective resistance increased and this A/ 
contrast was again easily observed. 

In order to investigate the tip-sample distance, the bias 
was reduced to zero partway through image acquisition. 
The result is shown in Fig.jsj^b). In this case, the contrast 
in A/ disappears with the lack of an applied bias voltage. 
Atomic contrast in A/ should be observable without an 
applied bias voltage [19]. In order to observe A/ con- 
trast with no applied bias voltage, we needed to advance 
the tip 340 pm closer to the surface. The remainder of 
the image, in Fig. |3|c) on, clearly shows the expected 
contrast in A/. 

One might ask what effect removal of the bias volt- 
age has upon the average tip-sample distance, knowing 
that there is an electrostatic force between tip and sample 
that scales as the square of the voltage difference [20] . We 
have found, in agreement with previous studies [21], that 
a model of the tip that incorporates a spherical apex with 
a conic structure to account for long-range electrostatic 
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FIG. 3. Top: Schematic of the experiment. Bottom: Simul- 
taneous / and A/ data acquired at 4.2 K. a) Vtip = —1.5 V. 
b) Vtip = OV. c) Tip is approached 340 pm closer to the sur- 
face. Data were collected at A = 100 pm, /o = 16 777 Hz and 
images are 8nmx8nm. 



interactions and a small nano-tip with negligable electro- 
static contributions but with which imaging is performed, 
as shown in Fig. [sj^a), is quite accurate when describing 
long-range electrostatic forces. Typical values to describe 
our tips are rtip = 5nm, a = 70°, and z^^^ = 1 nm, re- 
sulting in an average force, over one oscillation, of 30 nN 
when Vtip = 2V. Following p!8l, the effect of removing 
this bias is that the average tip-sample distance would in- 
crease by only 17 pm. Thus, even accounting for the effect 
of removing the bias voltage, we needed to approach the 
tip over 300 pm to the surface to observe the A/ con- 
trast in Fig. [sjc). It is therefore not possible that the 
frequency shift increase (e.g. in Fig. [l] or Fig. |3ja)) was 
recorded in region I where Pauli repulsion dominates. 

We now consider the effect of the sample having a re- 
sistance given by Rg^ with the tip biased at Vtip- The 
electrostatic interaction between tip and sample causes a 
force that scales as the square of the voltage difference, 
CK;ip — "^ampie)^- We ucglcct the coutact potential differ- 
ence, Kpd, because it is simply an offset of the applied 
bias voltage. While local variations in Vc^d have been 
reported on the 7x7 surface with Kelvin probe measure- 
ments [22 , they would not lead to an increase in A/ over 
adatoms in both bias polarities. Letting K represent the 
prefactor which is independent of bias [20 : 
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FIG. 4. Top: Schematic of the experiment. Middle: Simul- 
taneous / and A/ data. Bottom: A/ versus / data, a) Data 
was acquired with no additional resistor between sample and 
virtual ground, b) = 10 MQ was inserted between sample 
and ground, c) R2 = 30 Mil was inserted. See text for de- 
tails. Data were acquired with A — 400 pm, /o = 19 390 Hz 
and images are lOnmxTnm. 
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A/ is, in the small amplitude approximation, directly 
proportional to kts^ and if k^^ is dominant. A/ would 
have a linear response with respect to / and a slope that 
is linearly proportional to Va^. 

To test this theory, we performed an experiment shown 
schematically in the top third of Fig. |4] A switch was in- 
stalled between the high-doped sample and the virtual 
ground. The virtual ground is provided by the opera- 
tional amplifier used to detect the tunneling current. The 
switch is used to add known resistances = 10 M17 and 
R2 = 30 Ml]. A/ should now follow the relation: 



M = h^mV,,^{Rs^R)I 



(3) 



F^^ = —K (Fjp — 2 Ftip I Rs ^ RI) (1) where R is either 0, Ri or R2, depending on the switch, m 

is independent of Rg^ Vtip and /, and b is independent of 



Assuming ^sample = I Rs <^ Vtip, the second order term 
can be discarded. Proceeding with the derivative to k^g. 



/. Considering the /-dependent term in Eq. [3j m and Rs 
are unknown. The data collected ai R = and R = Ri 
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can be used to determine m and Rg , and then to predict 
a slope of -0.99 Hz/nA when R = R2. The data and the 
slopes for the three switch settings are shown in Fig. |4] 
the observed slope when R = R2 is -1.03 Hz/nA. 

This strong agreement with our simple model, incor- 
porating a sample resistance Rs^ prompts us to further 
explain the mechanism of this model. Our fit parame- 
ters, for example, indicate that for the low-doped sample 
Rs = 164 Ml]. While this number seems high, the tun- 
neling current is injected in a very small area, and that 
this in effect creates a large current density that must 
disperse through a relatively small area. Assuming that 
the current / is injected in an area of radius 1 A, and 
then disperses through the sample radially, we can use 
classical electrodynamics to determine the order of mag- 
nitude of the expected resistance. The current density 
then scales as 1/(2 7rr^), where r is the distance from 
the current injection, and is directly proportional to the 
derivative of the electrochemical potential via the sam- 
ple resistivity. The voltage between the area in which 
the tunnel current is being injected and a position in the 
bulk sufficiently far is thus 
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Given this sample has a resistivity p = 6 to 9Qcm, Rg 
is in the range of 96 to 143 Ml], which agrees well with 
the fitted Rs value. In this simple picture, the voltage 
drop would be highly local, but the charge density and 
electric field are both most intense near the tip apex, 
where one would expect the highest contribution to the 
electrostatic attraction. This picture also has a much 
higher error when we consider the high-doped sample, 
and it is likely that an atomic-scale theory of electronic 
conductance is required to fully explain the observed Rs . 

To summarize, we observe a frequency shift increase 
over adatoms that scales linearly with current. We have 
characterized this response to current as a function of 
bias voltage. This frequency shift can be explained by the 
effect that the tunnel current has on the surface potential, 
given non-negligable sample resistivity. 

As far back as simultaneous FM-AFM and STM have 
been attempted, contrast inversion has been observed as 
function of tip state and of applied bias voltage [THS]. 
Guggisberg and coworkers suggest that short-range elec- 
trostatics might explain the contrast inversion, but do 
not propose a model for this [10] . Contrast inversion can 
be explained within our model quite easily: At low bi- 
ases, the tip images the adatoms as predicted by theory, 
and shows a negative A/ contrast over adatoms, while 
at higher biases, the A/ contrast is due primarily to the 
decrease in capacitive attraction due to the sample resisi- 
tance and tunnel current, as reported in this Letter. 

High spatial resolution demands short-range forces, 
which implies small tip-sample distances. Given an inter- 



est in simultaneous AFM and STM, one must be aware 
that a moderate tip-sample bias can also produce a tun- 
nel current so large that it would affect the surface po- 
tential and measurable A/ contrast. For surfaces such as 
Cu, with a p ~ 2 X 10~^ l^m, this effect is ignorable [6j; 
even at Vtip = lOV and / = 100 nA, the expected A/ 
would be < 1 mHz. However this effect must be taken ac- 
count of when performing combined STM/ AFM of any 
surface with appreciable resistivity, including, of course, 
semiconductor surfaces. 

We thank J. Repp, M. Emmrich and M. Schneider- 
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